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Abstract 



It is shown that the quark mass density-dependent model can not be used 
to explain the process of the quark deconfinement phase transition because 
the quark confinement is permanent in this model. A quark mass density- 
and temperature-dependent model in which the quark confinement is imper- 
manent has been suggested. We argue that the vacuum energy density B is 

. T' 

a function of temperature and satisfies B = Bq 



1 



where T r is the 



critical temperature of quark deconfinement. The dynamical and thermody- 
namical properties of bulk strange quark matter for quark mass density- and 
temperature-dependent model are discussed. 

PACS number: 12.39.Ki,21.65.+f,25.75.Dw,11.10.Wx 
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I. INTRODUCTION 



It is generally accepted that the fundamental theory of strong interaction is quantum 
chromodynamics (QCD), however in reality, because of its nonperturbative characters, QCD 
has very little impact on the study of low and medium energy nuclear phenomena. Many 
effective models, some of them based on quark and gluon degrees of freedom and the 
others based on nucleons and mesons f2j, or quarks and mesons HU, have been employed 
to investigate the nuclear matter and strange matter. The quark mass density-dependent 
model (QMDD) suggested by Fowler, Raha and Weiner ||, is one of such candidates of 
effective models. Though this model involves some arbitrary choices and can not reproduce 
all conclusions given by lattice calculations, it is introduced as an alternative to the static 
bag model of confinement and has substaintial application in the study of bulk quark matter, 
especially strange quark matter and strange star. 

Recently, since the speculation of Witten || that the strange quark matter (SQM) may 
be more stable than normal nuclei, especially, since the argument given by Greiner and his 
co-workers that the small lumps of strange quark matter (strangelets) may be produced 
in relativistic heavy-ion collisions and could serve as an unambiguous signature for the 
formation of quark-gluon plasma, much theoretical effort has been devoted to studying the 
properties of SQM. Many investigations have been carried out in the frame of MIT bag model 
IHJ or QMDD model. Obviously, a successful effective model should be used to describe not 
only the dynamical and thermodynamical properties of SQM, but also the phase transitions 
of QCD. 

In this paper, we will focus our attention on the QMDD model. According to the QMDD 
model, the masses of u, d quarks and strange quarks (and the corresponding anti-quarks) 
are given by 



B 

m. 



q 3n B ' 



(q = u, d, u, d), (1) 



B 

m S! s = m s0 + - — , (2) 
3n B 

where n B is the baryon number density, m s0 is the current mass of the strange quark and 
B is the vacuum energy density inside the bag. At zero temperature, 

n B = \{n u + n d + n s ), (3) 

n u ,rid,n s represent the density of u-quark, d-quark and s-quark, respectively. The basic 
hypothesis Eqs.([J) and (fj) in QMDD model can easily be understood from the quark con- 
finement mechanism. A confinement potential which is proportional to r (or r 2 ) must be 
added to a quark system in the phenomenological effective models because the perturbative 
QCD can not give us the confinement solution of quarks. The confinement potential kr pre- 
vents the quark goes to infinite or to the very large regions. The large regions or the large 
volume means that the density is small. This mechanism of confinement can be mimicked 
through the requirement that the mass of an isolated quark becomes infinitely large so that 
the vacuum is unable to support it. Thus, for a system of quarks at zero temperature, the 
energy density tends to a constant value while the mass tends to infinity, as the volume 



2 



increases to infinity or the density decreases to zero fllQ| . This is just the picture given by 
Eqs.(jl])-(||). In fact, the similar confinement mechanism emerges in MIT bag model also. 
The boundary condition of confinement for MIT bag corresponds to that the quark mass is 
zero inside the bag but infinity at the boundary or outside the bag [11 



Although the QMDD model can provide a dynamical description of confinement and 
explain the stability and many other dynamical properties of SQM at zero temperature, 
when we extend this model to finite temperature and discuss the thermodynamical behaviors 
of SQM, many difficulties will emerge. Firstly, the thermodynamic potential Q is not only 
a function of temperature, volume and chemical potential, but also of density, because the 
quark masses depend on density. How to treat the thermodynamics with density-dependent 
particle masses self-consistently is a serious problem and has made many wrangles for this 



model in references [I2|-|T5| . Secondly, as will be shown below, it can not reproduce a correct 



lattice QCD phase diagram qualitatively or give us a successful equation of state when 
ris —> 0. It can not describe the phase transitions of quark deconfinement because the 
quark masses are independent of temperature. To overcome this difficulty, we will suggest a 
quark mass density- and temperature- dependent model (QMDTD) in this paper. Instead 
of a constant B in Eqs.(|l|) and (0), we argue that B would be a function of temperature and 
choose the function B(T) from Friedberg-Lee model. We will prove that the above difficulty 
can be overcome in our QMDTD model. 

The organization of this paper is as follows. In the next section we review three different 
treatments concerning the thermodynamics with density-dependent quark mass in refs. |12| 



T3| and |14| respectively, and prove all treatments can not overcome the second difficulty 



mentioned above. In section 3, we give detailed arguments on the temperature dependence 
of vacuum energy density B and extend QMDD model to a QMDTD model. Our results 
are summarized in section 4. In this section we prove that the temperature T vs density % 
phase diagram for QMDTD model becomes reasonable and it can be employed to mimic the 
QCD phase transition qualitatively. The comparison of QMDD model and QMDTD model 
for studying the dynamical and thermodynamical properties of SQM will also present in 
section 4. The last section is a summary. 



II. THERMODYNAMICAL TREATMENTS 

At finite temperature, the antiquarks must be considered. Eq.([|) becomes 

n B = i(An u + An d + An s ), (4) 

where 

An =n _ n = 9i r d 3 k ( i i 

i Ul Ul (2vr) 3 io \exp{/3(£i - m)] + 1 exp[/3( £l - + ^)\ + 1 

(nj) rii is the number density of the (anti)flavor i {i = u,d, s), gi = 6 is the degeneracy 
factor, Hi is the chemical potential (for antiparticle m = —fii). Inside SQM, s (and also s) 
quarks are produced through the weak processes 

u + d < — > u + s, s — > u + e~ + z/ e , d — > u + e~ + u e , u + e~ — > d + u e , (6) 
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and similarly for antiquarks. The system of SQM must satisfy the following constraints. 
The condition of chemical equilibrium yields ]12| 



Us = Vd, Hs = Hu + He- (7) 

The condition of charge neutrality reads 

2An u = An d + An s + 3An e . (8) 
The thermodynamic potential of SQM system is 

n = E n i = -E(^8 / o °°^ln(l + e-^-^), (9) 

where i stands for u,d,s (or u,d,s ) and the electron e(e + ), gi — 2 for e and e + . Not- 



ing that Si(k) = ymj + k 2 and m 8 is given by Eqs.(|T]),(@) and (|), we can calculate the 
thermodynamic potential Q under the constraints Eqs.(^,(M) and (^j), and obtain the ther- 
modynamical quantities such as number density rii, pressure, internal energy and etc.. Due 
to the density-dependent quark mass, many different treatments had been given in the 
references. 

A. First treatment 



The first thermodynamical treatment for QMDD model was given by Chakrabarty ||T2 . 
After getting the thermodynamic potential Q, he used the usual thermodynamical formula 
to calculate the number density rii, total pressure p and the total energy density e, and 
found 



1 dQ 

V djii 



(10) 

T,n R 



P=~y, (11) 



W -yQf 



(12) 



After comparison with the results given by MIT bag model, Chakrabarty claimed the prop- 
erties of SQM given by QMDD model were found to be very different from those predicted 
by the MIT bag model. Since the density dependence of quark mass has not completely 
and explicitly taken into account in this thermodynamical calculations, this treatment seems 
incorrect fl3|]. But in order to compare with other treatments, we list this treatment here 
also. 
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B. Second treatment 



The second different treatment is given by Benvenuto and Lugones ||13|| . They claimed 
that the features found by Chakrabarty are consequences of an incorrect thermodynamical 
treatment for QMDD model. In deriving the energy density and the pressure, an extra term 
appears due to the dependence of the quark mass on the baryon density. 

The results become 



p = - 



1 d(Q/n B ) 



V d(l/n B ) 



Q n B dVL 



(13) 



(14) 

and rii still satisfies Eq.(|10|). The extra term produces significant changes in the energy per 
baryon, make the pressure take the negative value in the low density region and shift the 
stability window of strange matter (SM). In almost all cases they found that the properties 
of SQM in the QMDD model are nearly the same as those obtained in the MIT bag model. 



C. Third treatment 



The third different treatment is done by Peng and his coworkers [H. Their improvements 
include: (1) Based upon a quark condensates argument, instead of Eqs.([l|) and @, they 
introduce 



D 

T/3 ■ 

B 



(q = u, d, u, d), 



n 



m. 



m s0 + 



D 

T/3- 
13 



(15) 
(16) 



n 



where D is a parameter usually determined by stability arguments. (2) They agree with the 
second treatment that one must add an extra term to the pressure formula because of the 
quark mass density-dependence, but do not agree with them for adding an extra term to 
the expression of the energy density because it can not give a correct QCD vacuum energy. 
The pressure and the energy density given by this treatment are 



P = ~77 



1 d(Q/n B ) 



V d{l/n B ) 



Q n B dVt 



(17) 



7>i 



SI ^ T dVL 



(18) 



In fact, this treatment is a "mixture" of the first and the second treatment. It chooses the 
pressure of the second treatment and the energy density of the first treatment as its pressure 
and energy density respectively. 
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Now we are in the position to study the thermodynamical behavior of SQM by us- 
ing the QMDD model. The temperature T vs density curves are shown in Fig.l 
by three dashed lines for three treatments respectively where we choose the parameters 
B = 170MeVfm" 3 , m s0 = 150MeV, D = HOMeVfrn" 1 and P = 400MeVfm" 3 . We see 
from Fig.l that the temperature T tends to infinite when % — > 0. This result is treatments- 
independent and can easily be understood if we notice the basic hypothesis of QMDD model, 
namely, the Eqs.([IJ) and (|^) (or Eqs.flTS|) and (P]|)), the quark masses are divergent when 
tib — > 0. To excite an infinite weight particle, one must prepare to pay the price of infinite 
energy, i.e. infinite temperature. This result demonstrates that the confinement in QMDD 
model is permanent. The quark can not be deconfined for any temperature. This model can 
not describe the quark deconfinement phase transition and give us a correct phase diagram 
of QCD. 



III. QUARK MASS DENSITY- AND TEMPERATURE- DEPENDENT MODEL 

Obviously, if we hope to employ the QMDD model to mimic the phase transition of 
QCD, the first problem is to avoid the permanent confinement mechanism given by Eqs.(^) 
and (E3) (or Eqs. (15) and (IE)). It would be useful to recall what happen in MIT bag model 



and Friedberg-Lee soliton bag model ||16|| . MIT bag model is a permanent quark confine- 
ment model because the confined boundary condition does not change with temperature. 
The vacuum energy density B is a constant in MIT bag model. Contrary, the Friedberg-Lee 
soliton bag model is a impermanent quark confinement model. Its confinement mechanism 
comes from the interaction between quarks and a non-topological scalar soliton field. Since 
the spontaneously breaking symmetry of scalar field will be restored at finite temperature, 
the non-topological soliton will disappear and the quark will deconfine at critical tempera- 
ture. In this model, the vacuum energy density B equals to the different value between the 
local false vacuum minimum and the absolute real vacuum minimum. This value depends 
on the temperature. It means that B must be a function of temperature in Eqs.(ffl) and (0) 
if we hope to deconfine quark. In order to describe the phase transition of QCD, we must 
extend the QMDD model to a quark mass density- and temperature- dependent (QMDTD) 
model and suppose that B is a function of temperature. This is our first argument. 

Our second argument comes from the calculations of the effective masses of nucleons 
and mesons recently. We can sum the tadpole diagrams and the exchange diagrams for 
mesons by Thermo Field Dynamics and find the masses of nucleons and mesons all decrease 
with temperature [|17] -|20||. This result for p-meson is in agreement with recent experiments 



PT| , |T9| ]. According to the constitute quark model, the nucleon are constructed by three 
quarks and the meson by two quarks. It means that in a satisfying quark model we must 
consider the temperature dependence of the quark mass. But this effect has not been taken 
into account by Eqs.(|T]) and (fj) if B is a constant. 



According to the conclusion of Benvenuto and Lugones [|T^], the results found by QMDD 
model are nearly the same as that obtained in the MIT bag model. On the other hand, as 
was pointed out by []16 |, the MIT bag model can be obtained from Friedberg-Lee soliton 



bag model provided fixed some parameters. Then it is natural to choose B{T) given by 
Friedberg-Lee model as our input. 
Introducing an ansatz 
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B = B 



, < T < T c 



B = 0, T > T, 



(19) 
(20) 



where T c is the critical temperature of deconfmement. When < T < T c , Eqs. ([I]) and (0) 
become 



2>n 



B 



rp x 2' 
1- ( — 



T, 



(g = u, d, u, d), 



B 

3n B 



1 - { T C 



(21) 
(22) 



The masses of quarks not only depend on the density % i but also on the temperature T. 
And when T > T c , m q = 0, m S)S = m s0 . When T = 0, Eqs.(pT]) and fl22"|) reduce to Eqs.(P 
and @, and our QMDTD model reduces to QMDD model. 

In our later calculations, we prefer to the thermodynamical treatment of Eqs. (|T7j) and 
(|18|), because they consider the density- dependent masses and the QCD vacuum energy ex- 
plicitly. Substituting Eqs.(|2T|) and (|22| ) into Eq.(^), under the constraints Eqs.(|^),(0) and 
(|), and using the same argument as that of the third treatment, we find the thermodynam- 
ical quantities rii, p, e, which is still expressed by Eqs. ([It]), ([H]) and ( |T8D because Eqs. fl2~T|) 
and (|22"D have the same density- dependence as that of Eqs.(jl]) and (|j). Even though the 
expressions of rii, p and e are the same, but we would like to emphasize that the results 
given by our model and QMDD model are different because the thermodynamical potentials 
calculated from Eqs.(^) and (g) and Eqs.(|T9"D-(|2"2"D are quiet different at finite temperature. 
Our results are summarized in next section. 



IV. RESULTS AND DISCUSSION 



The numerical calculations have been done by adopting the parameters B = 
170MeVfm -3 , m s0 = 150MeV, T c = 170MeV. The temperature T vs baryon number density 
tib figure is shown in Fig.l where the pressure P is fixed to be 400MeVfm -3 . The three 
dashed lines refer to three different treatments of QMDD models, respectively, and the solid 
line refers to our QMDTD model. We see from Fig.l that the basic difference between 
QMDTD model and QMDD model is: when % — > 0, the temperature approaches to a 
critical temperature T c = 170MeV in our model, and diverges in QMDD model. It means 
that QMDTD model is a impermanent confinement model. It can be employed to describe 
the phase transition of QCD qualitatively. 

The same curves of QMDTD model but for different pressures p = 400, 300, 200 and 
150MeVfm -3 are shown in Fig. 2 respectively. We see from Fig.2 that the conclusion T 
approaches to T c when n B — > will not change with pressure. The basic physical reason is 
that in QMDTD model, B is a monotonously decreasing function in the region < T < T c 
and becomes zero when T approaches to T c . The singularity of quark mass at zero density 
of QMDD model has been wiped out in QMDTD model. When T = T c , m 9i9 - = and 
m s,s = m so. It guarantees that the divergence of temperature at zero density will not happen 
in QMDTD model no matter the values of pressure be. 
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To compare our model with QMDD model further, we investigate the thermodynamical 
stability of SQM at finite temperature. The energy per baryon vs baryon density curves at 
T = 50MeV for QMDTD model and QMDD model are shown in Fig. 3 where the solid line 
refers to QMDTD model and the three dashed lines for three different treatments of QMDD 
model respectively. We see that the solid line is lower than the others. It means that the 
SQM described by QMDTD model is more stable than that by QMDD model. The values 
of ubo and (e/n^o at the saturation point are summarized at Table 1. We see from Table 1 
that the value of ubo for QMDTD model is situated between the maximum value 0.55fm -3 
and the minimum value 0.36fm~ 3 for the second and the third treatments of QMDD model, 
but the energy per baryon (e/r^o = 1006MeVfm~ 3 is lower than all treatments. The point 
marked with a heavy dot in the solid line is the zero pressure point for QMDTD model, as 
can be seen clearly, which matches the lowest-energy state and satisfies the basic requirement 
of thermodynamics pointed by ref. |L4||. The results obtained from Fig. 3 represent that our 
model can reproduce all thermodynamical properties of SQM which has been explained by 
QMDD model. 

The study of the equation of state for QMDTD model will show that this model is 
suitable for describing the thermodynamical behavior of SQM. The curves of pressure vs 
energy density are shown in Fig.4 where the solid line refers to QMDTD model and the 
dashed lines for different treatments of QMDD model respectively. We see from Fig.4 that 
the behavior of solid curve is very similar to that of the second treatment. It is monotonous. 
The values of pressure become negative at low density region and asymptotically tends to 
the ultrarelativistic case at high density as expected, because of the asymptotic freedom of 
quark [13| . 

Finally, we hope to investigate the so called "stability window" of SQM at zero temper- 
ature ||13|| . According to the argument of Farhi and Jaffe the conditions under which 
the strange matter be a true hadronic ground state read: at P = 0, E/ns < 930MeV for 
strange matter and E/ub > 930MeV for two flavor quark matter. Noting that even at 
zero temperature, the formulae of QMDTD model are still different from that of the sec- 
ond and the third treatments of QMDD model, because instead of Eqs.(|l5|), flilf ) for the 
third treatment, and Eq.(|14]) for the second treatment, we have Eqs. ([]]), (|2|) and Eq. (|18l), 
respectively. Our result is shown in Fig. 5, where for comparison, the stability window of 
the second treatment of QMDD model is also plotted. We see from Fig.5, the regions of 
Bq for stable quark matter are different for these two cases. Bo is limited narrowly in 
69.05MeVfm~ 3 < B < 111.6MeVfm~ 3 for the second treatment of QMDD model, but 
widely in 168.7MeVfm~ 3 < B < 273.3MeVfm -3 for QMDTD model. The stability window 
is still trianglelike but the adjusted parameters B and m S Q can take more values for which 
the system is stable in our model. 



V. SUMMARY 

In summary, it is found that the QMDD model can not be used to describe the quark 
deconfinement phase transition because the temperature diverges when baryon number den- 
sity approaches to zero. The quark confinement in this model is permanent. In order to 
overcome this difficulty we suggest a QMDTD model in which the quark confinement is 
impermanent. We argue that the vacuum energy density inside the bag B be a function 



S 



of temperature and prove that the divergence difficulty of temperature dose not emerge in 
QMDTD model. This model can mimic phase transition of SQM qualitatively. Finally, 
we compare the dynamical and thermodynamical properties of QMDTD model with three 
treatments of QMDD model, and find that our QMDTD model is useful to describe the 
properties of SQM. 

This work was supported in part by the NNSF of China and the Foundation of Education 
Ministry of China. 
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VI. TABLE 



Table 1, The value of saturation point for different models: 







(e/n B ) MeVfm- i 


1st treatment of QMDD model 


0.46 


1023 


2nd treatment of QMDD model 


0.55 


1083 


3rd treatment of QMDD model 


0.36 


1120 


QMDTD model 


0.45 


1006 



VII. FIGURE CAPTIONS 

Fig.l The temperature T as a function of baryon density n# with a fixed pressure P = 
400MeVfm~ 3 , three dashed lines are for the first, second and three treatments of QMDD 
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models respectively, and the solid line is for QMDTD model. 

Fig. 2 The temperature T as a function of baryon density % for QMDTD model with 
four different fixed pressure P = 400, 300, 200 and 150MeVfm -3 respectively. 

Fig. 3 The energy per baryon e/n B as a function of baryon density n B for QMDD model 
(dashed line), and for QMDTD model (solid line). 

Fig. 4 The pressure P as a function of energy density E/V for QMDD model (dashed 
line), and for QMDTD model (solid line). 

Fig. 5 The stability windows of SM for the second treament of QMDD model and for 
QMDTD model. 



11 







1st treatment 

of the QMDD model 



3rd treatment 
of the QMDD model 



treatment 
QMDD model 




n B (fm ) 



FIG.1 



P=400MeVfm 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 

n R (fm" 3 ) 



FIG.2 



900- 
800- 
700- 
600- 
500- 
400- 
300- 
200- 
100- 

o- 

-100- 



1 st treatment 
of QMDD model 




3rd treatment 
of QMDD model 



QMDTD model 



-200 



o 



500 



1000 



1500 



2000 



E/V(MeVfm ) 

Fig.4 



2500 



220- 
200- 
180- 



stability window 
for 2nd treatment 
of QMDD model 



stability window 
for QMDTD model 




80 100 120 140 160 180 200 220 240 260 280 

-3, 



B(MeVf m ) 

Fig.5 



